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Abstract. Using a multi domain spectral method, we investigate systematically the general-relativistic model for 
axisymmetric uniformly rotating, homogeneous fluid bodies generalizing the analytically known Maclaurin and 
Schwarzschild solutions. Apart from the curves associated with these solutions and a further curve of configurations 
that rotate at the mass shedding limit, two more curves are found to border the corresponding two parameter 
set of solutions. One of them is a Newtonian lens shaped sequence bifurcating from the Maclaurin spheroid 
sequence, while the other one corresponds to highly relativistic bodies with an infinite central pressure. The 
properties of the configuration for which both the gravitational and the baryonic masses, moreover angular velocity, 
angular momentum as well as polar red shift obtain their maximal values are discussed in detail. In particular, 
by comparison with the static Schwarzschild solution, we obtain an increase of 34.25% in the gravitational mass. 
Moreover, we provide exemplarily a discussion of angular velocity and gravitational mass on the entire solution 
class. 
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1. Introduction 

Various numerical methods have been developed to in- 
vestigate relativistic rotating models for extraordinarily 
compact astrophysical objects like neutron stars (Wilson 
119721 Bonazzola & Schneider 119741 Friedman et al. 

115551 Komatsu et al. I1989al I1989bl Lattimer et al. 
1T55U1 Neugebauer & Herold IT9"9"21 Herold & Neugebauer 
1T5521 Bonazzola et al. IT5531 IT55S1 Eriguchi et al. IT551 
Stergioulas & Friedman HMol Nozawa et al. IT55S}l . For 
reviews see Friedman 119981 and Stergioulas 1998 

The homogeneous and uniformly rotating star models 
were first studied by Butterworth & Ipser l|19751 11976(1 
who found that, in addition to the analytically known 
Schwarzschild and Maclaurin solutions, they are bounded 
by a sequence of configurations rotating at the mass shed- 
ding limit . The investigation of such limit configurations 
is instructive since certain physical quantities reach max- 
imal values there. Butterworth & Ipser for example esti- 
mated a 30% increase in mass for fixed central pressure 
owing to rotation. 

In the present Letter we completed their studies, find- 
ing two further limiting curves, in particular a sequence 
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1 Due to uniform rotation, a shedding of matter sets in when 
centrifugal forces balance gravity at the equator. Then a cusp 
at the star's equatorial rim appears. 



of stars with infinite central pressure and a sequence 
of Newtonian lens-shaped configurations that bifurcates 
from the Maclaurin spheroids before ending in a mass 
shedding limit (Bardeen 119711 Ansorg, Kleinwachter & 
Mcincl 2003b). All five limiting curves were found to cir- 
cumscribe entirely the general relativistic solution for ho- 
mogeneous star models that are continuously joined to the 
static Schwarzschild solution — hereafter called the 'gener- 
alized Schwarzschild solution'. This was done using the re- 
cently developed AKM method (Ansorg, Kleinwachter & 
Meinel 2002, 2003c), which allows one to solve the Einstein 
equations to high accuracy even for critical configurations. 
In particular we are able to determine to high precision 
the extreme configuration possessing both a mass shed 
and infinite central pressure. At this point several physi- 
cal quantities reach their global maxima and we provide 
explicit values for these. 

In what follows we use units in which the speed of light 
and Newton's constant of gravitation assume the value 1. 



2. Metric tensor and field equations 

The line element of a stationary, axisymmetric and asymp- 
totically flat spacetime describing a uniformly rotating 
perfect fluid body can be cast into the form 



ds 2 = e- 2U [e 2k (dp 2 + d( 2 ) + W 2 dcp 2 ] - e 2U (adcp + dtf 
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The corresponding Lcwis-Papapetrou coordinates 
(p, £, (p, t) are uniquely defined if we require continuity 
of the metric coefficients and their first derivatives at 
the body's surface, fn a comoving frame, for which the 
metric assumes the same form with metric functions U', 
k' , W' and a', the relativistic Euler equation can easily 
be integrated to determine the pressure p. For constant 
mass-energy density p this results in 



p = fi e 



y Q -u' 



1 



(1) 



where Vo is the constant surface value of U' . 

Taking into account asymptotic flatness, boundary and 
transition conditions at the surface and regularity along 
the rotational axis (p = 0), the interior and exterior 
field equations form a complete set of equations to be 
solved, which is done by applying the AKM method. For 
a comprehensive discussion of this multi domain spectral 
method see Ansorg et al. (2003c). 



3. Known static and Newtonian limits 

3.1. Schwarzschild solution 

Solving Einstein's equation for a static homogeneous star 
yields the famous Schwarzschild metric which reads in the 
above coordinates 



c = 



l-M/(2r) u 



l + M/(2r) 
for the exterior (r 



1 



MY 
2r) 



W = e k p a = 



1 



y/p 2 + C 2 > R) and 
- M/(2R) 1 - Mr 2 /(2R 3 ) 



l + M/(2R) l + Mr 2 /(2R 3 ) 



l + Mr 2 /(2R 3 ) H 

for the interior region (0 < r < R). Here M is the gravi- 
tational mass and R denotes the coordinate radius given 
implicitly through 



M = pf i? 3 



1 



MY 
2RJ 



For any physically relevant solution the pressure JIJ 
has to remain finite, which is fulfilled for R > M. This 
imposes an upper bound on the mass, namely 

4 1 



M < 



0.14477. . 



3.2. Maclaurin spheroids 

In Newtonian gravity the problem of self gravitating rotat- 
ing ideal fluids requires solving the Poisson equation for 
the body's gravitational field while satisfying the Euler- 
Lagrange equation governing its motion as an ideal fluid. 
This leads to a free boundary value problem. A particu- 
lar solution for uniform rotation and constant mass den- 
sity /i can be found by assuming the surface to be a 



spheroid. One obtains the so called Maclaurin spheroids, 
parametrized here by focal distance po and the ratio r p /r c 
between polar radius r p and equatorial radius r c . 

Having computed the gravitational field, the Euler- 
Lagrange equation is seen to be satisfied for a squared 
angular velocity 



il 2 = 2np£ [(3£ 2 + 1) arccot £ - 3£] , f = 



r 

1| - 1 



(bottom solid curve in Fig. 0. This relation holds inde- 
pendent of po- 

Note that Maclaurin spheroids exist for every r p /r 6 
[0,1] and po G [0, oo[, thus comprising a two parameter 
solution with arbitrary mass for fixed p,. 

3.3. First Newtonian lens sequence 

On the Maclaurin curve, an infinite series of points cor- 
responding to axisymmetric secular instabilities occurs, 
beginning at r p /r c = 0.17126, and accumulating in the 
Maclaurin disk limit r p /r c — * (Chandrasekhar I19fi7l 
Bardeen I1971|l . They are bifurcation points of further 
Newtonian sequences and correspond to singular post- 
Newtonian corrections (see Petroff l2003f) . The first one of 
these sequences is comprised of two segments that depart 
from the first bifurcation point. One segment proceeds to- 
wards the Dyson rings (Dyson IT£32l ITSMl Wong IWl 
Kowalewskv lTSgol Poincare ll^5alll885blll^5cl Eriguchi 
& Sugimoto 1981) whereas the other one ends in a mass 
shedding limit ( Bardeen IWT1 Ansorg et al. l2003b|l . The 
latter we will refer to as the 'first Newtonian lens sequence' 
motivated by the shape of the corresponding bodies. 

4. Generalized Schwarzschild solution 

For a systematic investigation of the general case the 
choice of parameters is not restricted by the numeri- 
cal approach. Nevertheless it is convenient to take non- 
ambiguous parameters that are restricted to a compact 
interval (here [0, 1]) in such a way that the endpoints rep- 
resent limiting configurations. Corresponding to the limits 
found, we selected the following magnitudes: 

— Mass shed parameter (as defined in Ansorg et al. 



d(p 2 



2 n-> 



r' p 



lim C 



dp 



where Cs(p) is the function describing the surface 
shape. Maclaurin and Schwarzschild bodies are char- 
acterized by (3 = 1 and the mass shed limit is fixed by 
(3 = 0. 

p c = p c /(p+p c ), where p c is the central pressure. Thus 
p c = stands for the Newtonian limit where mass 
and hence pressure vanish and p c — 1 for the limiting 
configurations with infinite central pressure. 
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Fig. 1. Squared angular velocity on the generalized 
Schwarzschild solution. Dotted lines indicate curves of 
constant central pressure (equally spaced in p c = p c /(M + 
p c ), p c being the central pressure and fi the constant mass 
density). Ergo-toroids appear for configurations situated 
above the dashed line. 
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Fig. 2. Magnification of Fig.Q]i n the vicinity of the first 
axisymmetric bifurcation point (C). Values for the first 
Newtonian lens sequence were taken from Ansorg et al. 
(2003b). Points (A) and (C) correspond to the cross sec- 
tions thus labelled in figure 6 ibid. 



Additionally we will use the ratio r p /r c of polar to equa- 
torial coordinate radius that is 1 only for Schwarzschild 
solutions. 

Our numerical analysis revealed that the generalized 
Schwarzschild solution is entirely bounded by the follow- 
ing five limiting curves (joined in the order listed): 

— the sequence of Schwarzschild solutions 
(r p /r e = 1, = 1, p c e [0,1]) 

— the Maclaurin sequence from the sphere to the first 
axisymmetric bifurcation point 

(r p /r c e [0.171, l],/3=l,p c = 0) 

— the first Newtonian lens sequence 
(r p /r c e [0.171,0.192], € [0,1], p c = 0) 

— a sequence of configurations rotating at the mass shed- 
ding limit (r p /r c e [0.192, 0.573], = 0, p c e [0, 1]) 

— a sequence of configurations with infinite central pres- 
sure (r p /r c e [0.573, 1], e [0, 1], p c = 1) 

This makes it possible to determine maximal values of all 
interesting physical quantities (see section 0. Moreover, 
we can now state that on the generalized Schwarzschild so- 
lution no quasistationary transition to a Kerr black hole is 
possible (in contrast to what was found for the relativistic 
Dyson rings by Ansorg et al. ( 2003a)) and that the surface 
remains convex in /^-coordinates. 

Fig. n depicts the (squared) angular velocity versus 
the radius ratio r p /r e . This is the completion of results 
of Butterworth & Ipser l|1975l IW7rJ|l . 2 It shows that the 
configuration with maximal angular velocity rotates at 
the mass shedding limit and possesses infinite central 
pressure. The magnification in Fig. demonstrates that 
the mass shedding curve does not terminate at the first 
axisymmetric bifurcation point (C), as was conjectured 

2 Note that in contrast to our work they used proper radial 
distances and kept the baryonic mass constant. 
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Fig. 3. Gravitational mass on the generalized 
Schwarzschild solution. The small interval on the 
r p /r e -axis, emphasized by a thick line, represents the first 
Newtonian lens sequence. Dotted lines indicate curves of 
constant central pressure (equally spaced in p c , see Fig. 

Ergo-toroids appear for configurations situated above 
the dashed line. 



by Butterworth and Ipser. Instead it is linked to the 
Maclaurin spheroids at this point via the first Newtonian 
lens sequence. 

The evolution of the gravitational mass for fixed cen- 
tral pressure is seen in Fig.[31and reveals a 34.25% increase 
in the maximal mass with respect to the static case. The 
global maximum is again found for the same configura- 
tion as for angular velocity. In both figures we included 
the line above which ergo-regions appear. Interestingly, 
this line corresponds roughly to one of constant mass. 

Other physical quantities as baryonic mass, angular 
momentum, polar red shift and circumferential radius 
show a similar behaviour: For fixed mass shed parameter 
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Physical quantity 


value 




Gravitational mass 


M = 0.19435 At" 1/2 


* 


Baryonic mass 


Mo = 0.27316 ^" 1/2 


* 


Angular velocity 


Q, = 1.8822 n 1/2 


* 


Angular momentum 


J= 0.03637^ 


* 


Polar radius 


r p = 0.04856 a* 1/2 




Equatorial radius 


r c = 0.08475 fj, 1/2 




Radius ratio 


r p /r c = 0.5730 




Circumferential radius 


R cilc = 0.41538 fi 1/2 


* 


Polar red shift 


Z p = 7.378 


* 



Table 1. Properties of the maximal mass configuration. 
Only valid digits are given. An asterisk indicates that the 
corresponding quantity has a global maximum there. 




Fig. 4. Meridional coss section (solid line) and ergo-region 
(dashed line) of the maximal mass configuration. Axes are 
scaled identically. 



or radius ratio they increase with increasing central pres- 
sure. Likewise they increase with decreasing mass shed 
parameter or decreasing radius ratio if the central pres- 
sure is kept constant. So a global maximum for each of 
them is obtained on the common edge of mass shed and 
infinite central pressure curves, corresponding thus to a 
very special limit star. 

5. Maximal mass configuration 

In Tabled masses and other quantities are listed for this 
extreme configuration. Because it resides on two critical 
curves we notice a loss in accuracy that was somewhat 
recovered by an extrapolation to an infinite approximation 
order m of the AKM method (cf. Ansorg et al. 2003c). 
Observe also the unexpectedly high value for the red shift 
Z p of a photon emitted at one of the poles. 

Fig. FJ] shows a meridional cross section of this con- 
figuration including the border of the ergo-region. As in 
general for diverging central pressure, the ergo-toroid de- 
generates by pinching together in the center. 

A detailed discussion, including the relation to rela- 
tivistic ring solutions, more realistic equations of state and 
further going aspects like stability will be published else- 
where. 
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